Longitudinal dynamics with RF phase modulation in the LNLS electron storage ring 
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In the Brazilian synchrotron light source, we observed that modulating the phase of the ac- 
celerating fields at approximately twice the synchrotron frequency suppressed remarkably well a 
longitudinal coupled-bunch mode of the beam driven by a higher order mode (HOM) in one of the 
radiofrequency (RF) cavities. In this work, we present the results of a set of systematic measure- 
ments, in single and multibunch mode, aimed at characterizing the effects of RF phase modulation 
on the beam. We compare those experiments with the results of tracking simulations and of a theo- 
retical model in which Landau damping is the stabilizing mechanism that explains the suppression 
of the longitudinal coupled bunch instability. We also measure the frequency of the stable islands 
created in longitudinal phase space by phase modulation and the longitudinal beam transfer func- 
tion (BTF) as a function of the modulation frequency and amplitude. The experimental results are 
in good agreement with theoretical expectations. 

PACS numbers: 29.20.Dh,29.20.-a,05.45.a 



I. INTRODUCTION 

Studies related to the effects of phase or amplitude 
modulation of the accelerating RF fields in circular ac- 
celerators date back to the early 1990's when a series of 
classical experiments were performed with the aim of un- 
derstanding basic aspects of non-linear longitudinal beam 
dynamics at low current in hadron machines 
Apart from the fundamental interest of those experi- 
ments, related to the possibility of direct observations of 
non-linear phenomena such as the formation of islands in 
phase space and measurement of synchrotron tune varia- 
tions with amplitude, the main applications envisaged by 
those works were related to issues such as lifetime lim- 
itations due to rf noise, the development of super-slow 
extraction techniques and the use of parametric feed- 
back for multibunch instabilities. Those early works laid 
the basic theoretical foundations for subsequent interest 
in applying RF phase or voltage modulation to electron 
storage rings used as synchrotron light sources 0, Q as 
a tool to control electron bunch density H, , increase 
beam lifetime and reduce the amplitude or even sup- 
press coupled-bunch instabilities [ToL ITTL IT^ . resulting 
in relevant improvement to the performance of those ma- 
chines. Some of those previous works analyzed the effects 
of RF phase modulation at the first and third harmonics 
of the synchrotron frequency 1 1(1 lla . , whereas others 
(such as the present paper) [llj focus on modulation at 
the second harmonic. 

However, in spite of the success of the implementation 
of such techniques at several laboratories, some aspects of 
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the physical mechanisms that explain the effectiveness of 
the reduction of the amplitude of coupled-bunch modes 
(CBMs) are still not completely understood and a quan- 
titative theoretical description is lacking. There is also 
some difficulty in simulating the effects of phase mod- 
ulation in multibunch mode; most of the work done so 
far simulates a single bunch and it is not obvious how to 
infer the multibunch behavior from these results. 

All experiments reported here were conducted at the 
LNLS 1.37 GeV electron storage ring ^5|, which rou- 
tinely operates to produce synchrotron radiation for 
users. The original motivation for our interest in RF 
phase modulation was the need to circumvent longitudi- 
nal instabilities that were observed in the LNLS ring right 
after the installation of a second RF cavity by the end of 
2003 0. This cavity has a longitudinal HOM (with a 
frequency of 903 MHz) which caused an intermittent or- 
bit distortion that mimics the second order dispersion 
function and which appeared when the beam coupled 
with the HOM of the cavity exciting a large longitudi- 
nal dipolar oscillation. With temperature and plunger 
scans we were able to identify the longitudinal mode LI 
(associated with the CBM JJ 133) of the new RF cavity 
as the main source of instabilities in the machine. Since 
it was not possible to find a cavity operational condition 
that would be free from instabilities, an active solution in 
the form of phase modulation of the RF fields at approx- 
imately twice the synchrotron frequency was attempted 
with success, having a noticeable impact on CBM ampli- 
tudes and helping to alleviate the orbit fluctuation. 

In this work, we develop a theoretical model and a sim- 
ulation tool to analyze the effects of phase modulation 
on the longitudinal beam dynamics. According to our 
model, RF phase modulation enhances stability through 
the increase of the incoherent frequency spread that re- 
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suits as the electrons populate the phase-space islands 
created by the action of the modulation. In order to ver- 
ify the model and simulation results, we perform a set 
of systematic measurements to characterize the effects 
of modulation on the beam dynamics observing directly 
the island tunes, the increased coherent damping and the 
changes to the longitudinal electron beam distribution 
caused by the modulation. 

The paper is organized as follows. We begin (Section 
II) with the motivation that led us to introduce phase 
modulation routinely in user shifts at the LNLS storage 
ring and with a report of our attempts to stabilize the 
beam after the installation of the new RF cavity. In 
Section III we calculate the electron equations of motion 
when subjected to phase modulation of the RF held at 
the second harmonic of the synchrotron frequency using 
the hamiltonian formalism and, using the Fokker-Planck 
equation, we derive the particle distribution in phase 
space. Also in Section III we calculate the dynamics for 
small oscillations around the stable fixed points created 
when modulation is turned on and, with the standard 
formalism used to calculate Landau damping, we find 
the increase in the amount of damping of centroid os- 
cillations related to phase modulation. In this section 
we also calculate the longitudinal beam transfer function 
and the stability diagram showing that the stable area is 
increased when phase modulation is turned on. In Sec- 
tion IV we present results of simulations of the effects of 
phase modulation in single bunch and multibunch modes. 
In those simulations one of the circulating bunches has 
an internal structure whereas all others are treated as 
macroparticles. Finally, in Section V, we compare the re- 
sults of theory and simulation with experiments in which 
we measure the island structure, the frequency for small 
oscillations around the stable fixed points created by RF 
phase modulation, the increase in the amount of damp- 
ing due to modulation and the longitudinal beam transfer 
function (BTF), verifying the results obtained with the 
theoretical model presented in Section III. Conclusions 
are given in Section VI. 



II. CBM STABILIZATION BY PHASE 
MODULATION 



TABLE I: LNLS ring main parameters. 
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shown in Figure ^ The distortion pattern has the same 
6-fold symmetry of the storage ring and shows up in the 
form of a fast transition between two states. For most of 
the user run periods, small orbit fluctuations with am- 
plitudes of ±3 /jm were observed and could be detected 
at the most sensitive beamlines. However, larger ampli- 
tude distortions were sometimes observed as shown in 
Figure ^ 

The origin of the orbit fluctuations was identified as 
being related to large amplitude dipolar longitudinal os- 
cillations caused by the interaction of the beam with 
HOMs excited in the cavities. In fact, the second order 
term of the non-linear dispersion function 771 (s) was mea- 
sured and has the same profile observed in the beam posi- 
tion monitors (BPM's) reading and shown in Figure ^3. 
When a longitudinal dipolar oscillation is present due to 
the influence of a HOM, the relative energy deviation is 
given by 5(t) = So cos il s t. The closed orbit distortion is 
then 

x e (s,t) = S(t) V (s) + 5 2 (t)r h (s). (1) 

Since the oscillation is fast (the synchrotron frequency 
£l s is of the order of 25 kHz) the BPMs measure only 
average positions 



The LNLS synchrotron light source is based on a 
1.37 GeV electron storage ring and operates with 250 mA 
of initial current in normal users shifts. A summary of 
the main parameters for the Brazilian machine is shown 
in Table HI At the end of 2003 the RF system of the ring 
was upgraded and a new cavity was installed |l6l Il7| . 
This second RF cavity was added to the system in order 
to prepare it for power demands that the future installa- 
tion of insertion devices will bring about. After operation 
with two cavities was initiated, the beam suffered from 
frequent instability outbreaks during user runs. 

The typical signature of the beam instability is an or- 
bit fluctuation measured only in the horizontal plane, as 



(x e ) = - Vl (s)6l (2) 

whereas the maximum amplitude of the oscillations is 
given by 



xr x ( S )=5ov(s) = \^iv(s). (3) 

x max snows U p m detectors of the beam lines as an 
effective increase of the horizontal beam size. The same 
orbit distortion can be obtained by phase modulating (at 
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Figure A 
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FIG. 1: Orbit distortion caused by a longitudinal CBM. The 
distortion occurs only in the horizontal orbit and has the same 
6-fold symmetry of the machine. (A) Time evolution of the 
horizontal orbit measured by four different BPMs during a 
user shift showing the orbit fluctuation caused by the insta- 
bility excited by a HOM. (B) Profile of the orbit fluctuation 
along the storage ring. This is obtained by subtracting the 
orbit right before the fluctuation from the orbit right after it. 



mode (CBM) instabilities as a function of temperature 
and plunger positions. By mapping the variation of the 
modes with these parameters, we determine the regions 
for which the growth rates of the oscillations are below 
the radiation damping rate, which for the LNLS storage 
ring at nominal operating energy is l/r rac ; w 270 s _1 , as 
shown in Figure [21 

The results of the HOM survey showed a comfortable 
situation for the old cavity. The temperature range for 
which there are no dangerous longitudinal HOMs is quite 
wide. Unfortunately, the same is not true for the new 
cavity. A list of some of its main longitudinal modes is 
shown in Table [H] Tc is the critical temperature of the 
mode, defined as the temperature at which the frequency 
of the HOM matches the frequency of the CBM oscilla- 
tion. The values of Ql have been measured in the new 
cavity but those of R s /Qq are simulated values. The 
main problem with the new cavity is the longitudinal 
mode LI whose critical temperature is just at the edge 
of the available temperature control range but is rather 
insensitive to changes of the cavity tuning parameters, as 
shown in Figure |21 An attempt was made to shift that 
mode using the plunger already installed in the cavity 
but the extent of the plunger motion was not enough to 
accomplish the task 
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TABLE II: Main TM 0m „ Modes of Cavity 2 
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the synchrotron frequency) the RF master generator that 
drives the RF stations. 

Attempts to eliminate the instabilities by changing the 
parameters of both cavities have proved to be not very 
effective. The search for a convenient operational point 
of the cavity was performed by scanning the parameters 
that can be used to control the cavity tuning (tempera- 
ture, plunger position and axial deformation) and mea- 
suring the corresponding HOM frequencies. 

The identification of the potentially harmful modes 
was possible by comparing them with the modes mea- 
sured for a similar cavity used in our injector booster 
synchrotron, which had been fully characterized in the 
laboratory, and by comparing the frequency shifts of the 
modes with those predicted from their simulated field 
maps. The procedure adopted to analyze the effect of 
the HOMs is the same used for similar cavities operat- 
ing at ELETTRA and ANKA [3 and consisted in the 
determination of the growth rates of the coupled bunch 
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FIG. 2: CBM growth rates for the new RF cavity in the 
best operation point for the plunger position, calculated for 
Eq — 1.37 GeV and / = 250 mA. Notice the absence of a 
mode free region. 

The coupled bunch mode associated with the LI mode, 
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CBM y 133 , is always present in the beam spectrum ex- 
cept in those regions where the growth rate of other 
modes, such as CBM ((126 associated to mode L3, are 
larger. HOM mapping was useful for the identification 
of the best operation point but this point turned out to 
be not good enough. As a result, an active solution in 
the form of phase modulation of the RF fields at approx- 
imately twice the synchrotron frequency was attempted 
with remarkable success. Figure [21 shows the longitudi- 
nal spectrum of the stored beam as measured by a button 
pick-up. The suppression of the dipolar synchrotron side- 
band of the harmonic #281 of the revolution frequency 
by more than 40 dB illustrates the damping of this CBM 
by RF phase modulation. 



FIG. 3: Measured frequency spectrum of the signal from a 
button pick-up in the storage ring demonstrating the sup- 
pression of dipolar synchrotron sidebands as a result of RF 
phase modulation at the second harmonic of the synchrotron 
frequency. 



III. THEORETICAL MODEL 
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A. Longitudinal Hamiltonian with Phase 
Modulation 



We initially follow the standard hamiltonian anal ysis 
for the longitudinal dynamics with phase modulation |l4| 
in which we expand the complete Hamiltonian near the 
second harmonic resonance of the synchrotron frequency 
and retain only the lower order terms. The Hamiltonian 
for an electron in a bunch subjected to phase modulation 
can be written as follows 



H(S,cf>;t) 



+ oj s tan (j)[sm (j) cos(A m sinui m t) + cos <p sm(A m sinuj m t)] 
u s cos 4> cos(A m sinuj m t) +u> s sm(j)sui(A m sinLLi m t) — Lu s t&n(j) s 

I 



(4) 



where uj s — 2nf s is the synchrotron frequency 
4> s = 7r — <j) s where <j) s is the synchronous phase, A m is 
the modulation amplitude and uj m is the modulation fre- 
quency. Using the action-angle variables (J, "J) we make 
the following transformation 



(5a) 
(5b) 



-V2Jsin(V> + to m t/2 + tt/4), 
^2J cos(-0 + u m t/2 + tt/4), 



which describes the trajectories of the particles in the 
phase space near the second harmonic resonance in a 
reference frame which rotates at half of the modulation 
frequency The resulting time averaged Hamiltonian K 
[Appendix |A"| is then: 



(K)t 



J 



16 



cos2?/> (6) 



where e = A rn tan <f) s 

The fixed points of the system are the solutions of 



dJ 

~dt 



dK 

dtp 



and 



~dt 



dK 
~dJ 







(7) 




Lu m < (2 + e/2)uj s 
uj m > (2 + e/2)oj s 



(8) 



corresponding to the angles ipSFP = 0,7r, which are sta- 
ble since cos 2-0 > 0, and the unstable fixes points are 



Jufp 



1 - 







- 2e , oj m < (2 - e/2)uj s 
, u m > (2 - e/2)u s 



(9) 



corresponding to the angles ipUFP = tt/2, 37r/2. 

The phase space coordinates of the new stable fixed 
points are related to the square root of Jsfp or Jufp 
and are functions of the RF phase modulation amplitude 
and frequency. 

The effect of phase modulation is to create new regions 
of stability inside the bunch, besides the original one. As 
the formation of those islands depends on the modulation 
amplitude and frequency, there is a continuum of forms 
from one island to the appearance of three islands when 
w m < (2 — e/2)u s and the formation of only two islands 
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FIG. 4: Phase space of the modulated bunch with two differ- 
ent modulation frequencies in the rotating reference frame, 
(a) uj m < (2 — e/2)io s with 3 stable fixed points and (b) 
uj m > (2 — e/2)ui s with 2 stable fixed points. The red line 
corresponds to the separatrix. 



when u> m > (2 — e/2)u> s . The longitudinal phase space 
for the situations described above is shown in Figure 01 



B. The Effect of Radiation 

All the arguments given above do not take into account 
the fact that electrons in a bunch lose energy through the 
emission of photons. This effect causes oscillations inside 
the bunch to damp, so it is important to know when this 
damping can destroy the formation of islands produced 
by phase modulation or modify some characteristics de- 
rived using the hamiltonian formalism. The analysis that 
follows uses a simple non-linear oscillation theory to cal- 
culate the trajectories of a single particle under the influ- 
ence of phase modulation at the second harmonic of the 
synchrotron frequency. 

The equation of motion for a particle subjected to 



phase modulation can be written as follows 



;ubjc 

m 



with 7^ the radiation damping time and 

1 



e = Am tan c 



and P=^l- 



(10) 



(11) 



We also consider that A m , the phase modulation ampli- 
tude, and d> are such that A m < cl < 1 rad and retain 



only the terms most important which are the parametric 
resonance term on the left side and the third order term 
on the right side of equation HlOjl . Non-resonant terms 
in A m and tan0 s are discarded. 

From the standard analysis of parametric resonances 
[2fj| we expect that the solution of equation I|1U|I has the 
form 



0(t) = a(t) cos (^t) + b(t) sin (^i) 



(12) 



Substituting 1121 in equation 1101 we can write a set of 
equations for the functions a and b 

2a + 2j d a + (a + 5u)b + kb{a 2 + b 2 ) = 0, (13) 



2b + 2 ld b + (a - 5u)b - ka(a 2 + b 2 ) = 0, 



where 



- ffl b - — 
a ~ 2 ' 8 ' 



Slu = 



- 2uj s 



(14) 



(15) 



If a particle is to be captured by the parametric reso- 
nance, the slowly varying coeficientes a(t) and b(t) must 
increase exponentially. Taking only the lower terms in a 
and b in equations l|13|) and (|14fl and using this condition 
we find that, in order for phase modulation to create sta- 
ble islands in the longitudinal phase space, it is necessary 
that the condition a > 2^ d be satisfied. Rewriting this 
condition in terms of the phase modulation amplitude we 
find that 1 



A™ > 



47d 



tan < 



0.026 rad 



(16) 



and this is a limit imposed by radiation damping bellow 
which no island is observed. 

Changing the longitudinal phase space reference frame 
to one which rotates with a frequency of oj m /2 in time we 
find that the coordinates in this new frame are related 
only to the functions a(t) and b(t) and the fixed points 
are given by the solution of equations ljT3|l and ifHJl when 
a = b = 0. Rewriting a — \/2rcos9 and b — \/2rsin(9 we 
end up with the following solutions 



-6U1 ± yjo? - 47^ 



2k 



tan 6 



2 ld 



(17) 
(18) 



If we set 7d = we find that the stable fixed points 
are the ones already calculated in the previous section 
using the hamiltonian formalism. To take into account 
the effect of radiation damping we define, using (|17fl . an 
effective amplitude of modulation which is modified by 
the radiation damping effect as: 



Aeff 
SI™ 



'A 2 



( 4 7d 



\ \ tan4> s \u> s 



(19) 



1 Calculated for a total accelerating voltage of 500 kV. 
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This quantity is particulary important when one uses the 
hamiltonian formalism to make predictions about mea- 
surements. In these calculations A e J n ^ should be used 
instead of A m when the damping effects are not negligi- 
ble. 



Dynamics of Small Oscillations Around the 
Fixed Points 



We now consider the dynamics of particle oscillations 
in the vicinity of the stable fixed points (SFPs) in or- 
der to derive two important properties of the electron 
motion: the equilibrium density distribution close to 
the SFPs and the frequency of synchrotron oscillations 
around those points (also known as the island tune), in- 
cluding its dependence on the oscillation amplitude. This 
is an important quantity since it describes the frequency 
of incoherent oscillation of the particles inside the island 
when subjected to a longitudinal kick and its dependence 
on oscillation amplitude is intimately related to the ef- 
fectiveness of phase modulation as a damping mechanism 
for coherent synchrotron oscillations. 

We start by expanding equation JSJ) around a fixed 
point with coordinates (Jo,?/>o) m the rotating reference 
frame [f| 



5' = — v2Jsin^ + y / 2Josin^?0) 
(/)' — V2J cos -0 — \J 2 J cos ipo, 



and find that 



H' = ^5' 2 + ^' 2 + 0(^) 



(20a) 
(20b) 



(21) 



where 



A = UJ S 



W s Jo UJ s Jq . 2 
a a Sm VO 



T> "III 



/ Jq Jo „ , , 

-— (1-- + T cos20 o ), 



W s Jo ^sJo 2 / 
a a C0S V0 



2 

UJs€_ 



■, Jo . Jo „ , 

1 + — + — cos 2ip 



(22) 



(23) 



and the rms bunch energy and phase spread are 

T 



as = 
where 

Id 



K 

Id 

1 

7~rad 



and 



250 s" 




and 



4 x lO -4 * -1 . 



(25) 



(26) 



In other words, particles close to the stable fixed points 
are subjected to an effective harmonic potential and we 
estimate the steady-state phase space profile of the elec- 
tron beam under phase modulation (when radiation ex- 
citation and damping are balanced) by a bi-gaussian dis- 
tribution. 

If we go further in expanding © we find that 



A -a , B ,12 UJ S ^2J S FP , /3 



16 



16 
24 



(27) 



and that the unperturbed frequency, or island frequency, 
is given by lj — V AB. This quantity is real for the stable 
fixed points and imaginary for the unstable ones. 
Changing variables to 




cos \fr' and 5' 




sin*' 



(28) 

and using the canonical perturbation technique |2l| with 
the generating function 

F 2 (*',/) = *7 + G 3 (/)sin31'' + Gi(/)sin* / , (29) 

where the functions G\ (I) and G3 (J) are chosen so that 
the term in J 3 / 2 is cancelled, we find that the phase 
averaged Hamiltonian is 



3 A 2 



B 2 e 

+ 



B 2 



8 A 2 B 2 

15 UJ 2 y/I S FP ( A 



4 A 2 B 2 

3/2 



1G 



B 



I 2 (30) 



and the revolution frequency around each stable fixed 
point for a particle with a phase amplitude <j> with respect 
to that point is 



In order to find the steady-state phase space profile of 
the electron beam under phase modulation around a sta- 
ble fixed point, we must solve the Fokker-Planck equa- 
tion, using H' from equation 121|) and taking into ac- 
count the radiation damping (7^) and quantum excita- 
tion (k = <r 2 jd) terms. When all forces are balanced the 
solution of the Fokker-Planck equation is a bi-gaussian 
where 



bj{4>) ~ lo 1 



3w s A 2 
~16 



B 



A 2 \B\ 



(31) 



where the parameters of the RF phase modulation are in 
the coefficients A and B. 



D. Beam Transfer Function and Landau Damping 



1 



27TCT5CT0 



exp 



1 

"to? 



1 



(24) To ex pl am why phase modulation has such an intense 

effect on the amplitude of CBM oscillations, we consider 
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that phase modulation creates a spread in the frequency 
distribution inside the bunch thus increasing the amount 
of Landau damping [23 • In this case we are only inter- 
ested in the damping of synchrotron oscillations of the 
bunch ccntroid which are produced by the coupling be- 
tween the synchrotron motion and a longitudinal HOM 
of the RF cavity. We therefore proceed to calculate the 
response of the beam centroid motion to an external har- 
monic perturbation when phase modulation is on. 

Suppose the particles in a bunch have a distribution 
^(r) in phase space. In case the beam is phase modu- 
lated the distribution can be approximated, in the rotat- 
ing reference frame, by three gaussian functions each one 
centered at a fixed point 



Mr) = 



Na 



2lTCT c 



N l2 

27T(7, 



(32) 



where the indexes c and il, 2 indicate the central island 
or the side islands formed when the phase modulation 
is turned on, r is a radial coordinate in the longitudinal 
phase space and Nj is the relative number of particles in 
each island such that Nn + Ni% + N c = 1. 

Considering that, when an external perturbation with 
amplitude Fq acts on the beam, a dipolar motion of the 
centroid is excited, we can write the electron distribution 
as 



(33) 



where V^i is a small perturbation. Using l|33|l and consid- 
ering that the equation of motion for each particle inside 
the island and near the SFPs can be approximated by 



t + 2 ld r + i? ci r = F e jnt 



(34) 



where r is the relative time displacement and v c .\ the 
island frequency, we can show that |2^| 



*i(r) 



N C F 



NgFg dt> 0tl 

uj c (fl — uj c ) dr uii(fl — u>i%) dr 
N i2 F dV 0t2 
lu c (Q — 0Ji2) dr 



(35) 



in which we introduced the damping rate due to radiation 
defining a complex frequency for each island: lo c ^ — v c ^ — 

iid- 

The centroid motion is given by the first moment of 
the phase space distribution 



f(fi) 



27T />00 





2uj c 



r 2 cos 6 'd6dr^(r, 9, t) 



N C I C {Q.) + N^Iiifl) 



(36) 



where the contribution of the side islands result in the 
same integral since their frequencies are the same. 




Frequency (kHz) 



FIG. 5: Calculated amplitude and phase of the beam transfer 
function for a bunch with (full red line) and without (dashed 
blue line) RF phase modulation and with half of the particles 
in the central stable island. 



The function I n (Cl) is the dispersion integral and is 
defined as 



I n (Q) 



r 2 dr d^o 
Q — uo n (r) dr 



(37) 



and the beam transfer function (BTF) for a beam under 
phase modulation in the rotating reference frame is 



J(fi) = N c I e ((l) + Ni-lli(Sl) 



(38) 



As the measurements are performed in the laboratory 
reference frame it is necessary to transpose equation (|38|) 
to this reference frame. To do this, note that in the 
lab frame the frequency of the lateral islands appears 
as a lower sideband of half of the modulation frequency 
[Appendix so that we have two different frequency 
responses: uj m /2+uj c for the central island and Lu m /2 — LUi 
for the side islands, as shown in Figure[3] It is interesting 
to observe that the effect of phase modulation is to create 
a second peak, related to the island frequency and also 
an intermediary phase jump in the BTF, these features 
are very different from the conventional picture expected 
when measuring the BTF of an electron beam without 
phase modulation. 

It is possible to estimate the effect of phase modula- 
tion on the damping time of the centroid motion if we 
calculate the increase in frequency spread due to modu- 
lation. The particle distribution in the frequency domain 
is given by 



N c 
Au; r 



+ 



Ni 



<<•>*-<•>)/ ^*H{ui - w), (39) 
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FIG. 6: Calculated distribution of incoherent frequencies in 
the rotating reference frame. (blue curve) without phase mod- 
ulation and (red curve) with phase modulation turned on, 
fm = 50 kHz, A m = 50 mrad and half of the particles in the 
central island. 



and shown in Figure El where H(x) is a step function. 
Defining the frequency spread as 



Auj = \/(uj 2 ) = U! 2 ^' q (lu) duj 



(40) 



we have that the ratio between the natural damping time 
to the one created when the beam is under phase modu- 
lation is 



7~nat 



Au nat /2ir + l/T rad 
/2w+l/T rad 



(41) 



The natural spread in the incoherent frequency 
for the bunch in the LNLS electron storage ring is 
(Auu)™^ w 160 Hz and the total damping, including the 
radiation damping time, related to this frequency disper- 
sion is T nat ~ 3 ms. Suppose half of the particles are in 
the central island, so that N c = 1/2, and the amplitude 
and frequency of the modulation are 50 mrad and 50 kHz 
respectively. The frequency spread and the new damping 
time due to phase modulation, given by equations l|40|) 
and (|TT|) . are Aui mo d = 490 Hz and r moi j = 2.2 ms. This 
damping time is smaller than the estimated growth time 
for the LI instability, which is 3.3 ms 2 . So, using the 
phase modulation we are able to set the beam well below 
the stability boundary and effectively damp this coupled 
bunch mode. 

So far we have only estimated the effects of phase mod- 
ulation on beam stability. We are going now to calculate 
the implications of the increase in frequency spread due 
to modulation on the beam when it is subjected to its 
own fields. 



2 This value for the mode growth time is for a temperature of 
40 °C, an accelerating voltage of 500 kV and a total stored cur- 
rent of 250 mA. 



The spread in incoherent frequencies in the bunch 
causes it to strongly damp coherent oscillations excited 
by a harmonic driving force. This force can be an ex- 
ternal excitation or wake field which interacts with the 
electron beam 22] . If we consider a short bunch which 
interacts with an impedance of the machine, such as a 
HOM, the centroid motion is governed by the following 
equation 



ds 2 



2 NrjUj r 



P 2 E a C 2 



(42) 



where s = vt with v = (5c is the particle velocity, N 
is the number of electrons in the bunch, e is electron 
charge, Lo rev is the revolution frequency, Eq is the parti- 
cles energy, C is the orbit length and Z\\ is the effective 
longitudinal impedance |22| . 
Introducing the ansatz 

f(s) =Be jas ' v , (43) 

the coherent synchrotron tune shift can be defined as 

je 2 Nr]LOrevC 2 



(ALu s ) coh = Q — uo 



2uj E q C 2 



(44) 



and the stability condition is so that: Im(Q) < 0. 

Therefore, using equations l|36|l and l|44|) we can define 
the stability diagram for coherent centroid oscillations in 
the following way 



U + iV = -j{Auj s ) coh = 



m g 2 (n) + p(n) 



(45) 



where /(w) and g(u>) are the real and imaginary parts of 
the BTF I (to) (equation 03). 

FigurcQshows the stability diagram for a bunch under 
phase modulation (black full line) and in normal condi- 
tions (dashed line). Note that, when phase modulation 
is turned on, the stable area is increased and that a new 
region, inside the loop, is formed. This new region is also 
unstable, as shown by the red line which is calculated 
for positive imaginary coherent frequency shift, what vi- 
olates the stability criterion. The dots shown in Figure[7| 
represent various operation points depending on the op- 
erating temperature of the new cavity. The green mark, 
for example, indicates the point for 50 °C that is unstable 
under normal operation conditions and becomes stable if 
phase modulation is turned on. 



IV. SIMULATION 

In order to gain some insight into the behavior of the 
bunches under phase modulation and verify the prop- 
erties derived from the theoretical model outlined above, 
we developed a longitudinal dynamics simulation code in- 
cluding beam loading effects and the presence of the lon- 
gitudinal HOMs as well as phase modulation. The track- 
ing code calculates the trajectories of individual particles 
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FIG. 7: Stability diagram in the U x V plane. (black full line) 
phase modulation on and for 50 % of electrons in the inner- 
most stable island, (blue full line) unstable line for a positive 
imaginary coherent frequency shift, (dashed line) phase mod- 
ulation off and (dots) various operating points for tempera- 
tures of the new cavity varying form 40 to 55 °C. 



A. Equations of Motion 

The model used in the simulation considers that ev- 
ery possible bunch is filled with electrons and each one 
of those bunches, except one, is represented by a single 
macroparticle. The bunches pass periodically through 
the accelerating cavities and contribute to the cumula- 
tive build-up of the induced fields. The effect of the beam 
loading is also taken into account in the calculation of the 
trajectories of the individual electrons of the structured 
bunch. 

In order to calculate the trajectories of these macropar- 
ticles in the longitudinal phase space, we use the following 
equations: 

U,n = Tl,n-1 ~ a^n-l^O, (46) 

S b , n = (1- 2 7d T )V-i + {eVtot - Uo \ (47) 

with 



in longitudinal phase space and derives the corresponding 
averaged quantities representing the motion of the bunch 
ccntroid as well as a projected (time averaged) longitu- 
dinal charge density profile appropriate for comparison 
with experimental results. 

Even though the application of the code to the single 
bunch case is relatively straightforward, the calculations 
for the multibunch case are much more time consuming, 
since at the LNLS storage ring we routinely operate with 
all 148 bunches filled (no gaps in the bunch train). In or- 
der to circumvent that difficulty, we use a hybrid simula- 
tion, in which we treat the first 147 bunches (&i — 6147) as 
macro-particles (no internal structure) whereas the last 
bunch (bo) is composed of 1000 individual macroparti- 
cles. Note that although each one of the 147 bunches 
without internal structure is represented basically by a 
single point in phase space (rather than by a distribu- 
tion), an effective bunch length (representing the projec- 
tion of the phase space bunch distribution on the time 
axis) is also attributed to those bunches when calculating 
the excitation of the wake fields. In fact, as the particles 
inside bo are redistributed in phase space as a result of 
phase modulation and interactions with the environment, 
its bunch length is also changed and that value of bunch 
length is then transferred to all non-structured bunches 
(b\ — 6147). In this way one important aspect of the inter- 
nal bunch structure is shared by all bunches. Clearly, this 
procedure cannot take into account all properties related 
to the internal motion of the bunches and describe their 
effect on the stabilization of the motion. However, we 
shall see in the following sections, that such simulations 
give additional indication that the increase in incoherent 
frequency spread inside the bunches is indeed the main 
stabilizing mechanism provided by phase modulation. 



Vtat = V rf +V b l (48) 

where r and S are the time displacement and energy de- 
viation of the macroparticle related to the reference syn- 
chronous particle, b is the index of the bunch and n is the 
number of turns, To is the revolution period, V r / is the 
gap voltage and V&z is the energy drained by the beam 
loading effect, which takes into account the beam loading 
of the fundamental mode and of the longitudinal mode 
(LI). In order to simulate each particle in the structured 
bunch (bo), we use the same equations and also consider 
radiation excitation. 

The phase modulation is included in the calculation of 
the rf voltage as follows 

V r f = Vo sm(<f> s - w r fn,n + A m sm(w m t)) (49) 

were Vo is the peak voltage, <p s the synchronous phase, 
Lu r f the angular RF frequency, A m the modulation am- 
plitude and oj m the modulation frequency which we set 
near the synchrotron second harmonic (uj m = 2u) s + 5u). 

The beam loading effect is introduced in the tracking 
simulation using the following expressions: 

V n = V n _ 1 e- u '^' M / 2QL+ju ^ M ~ kqoe-^r (50) 

where to res is the angular frequency and k is the loss 
factor of the mode, go is the charge per bunch and oy is 
the rms bunch length. 

B. Single Bunch Simulations 

Using the simulation code described above, we calcu- 
late the phase space distribution of the particles in a 
single bunch of 1000 macroparticles under phase modu- 
lation. Figure 03 shows the phase space with two different 
modulation frequencies and also without phase modu- 
lation. The phase space patterns observed are in good 



10 




phase modulation OFF - 
phase modulation ON 



FIG. 8: Single bunch (1000 macroparticles) simulation show- Frequency [Hz] 

ing the phase space for three different situations (blue) with- 
out modulation, (red) with f m < (2 — e/2)/ s and (green) FIG. 10: Amplitude of the dipolar oscillation of the bunch 
/ m > (2 — e/2)/ s . In the two last figures it is possible to centroid excited by a RF cavity longitudinal HOM for a single 
observe the formation of three and two islands respectively, bunch with 30 mA. 
as predicted by the theory. 




x [xNX ,z s] 



FIG. 9: Longitudinal profile in time domain of a simulation 
of 1000 macroparticles in three different situations: (blue) 
without modulation, (red) with f m < (2 — e/2)/ s and (green) 
f m > (2— e/2)/ s . Those patterns simulates the measurements 
of the signal of a button pickup. 



agreement with the ones calculated using the Hamilto- 
nian © and shown in Figure In particular, the posi- 
tion of the islands center given by simulation corresponds 
to the action amplitudes Js_fp=0.122 and 0.047 for 3 and 
2 islands respectively, as shown in Figure |H1 This is in 
very good agreement with the theoretical values 0.042 
and 0.130 for the cases with 3 and 2 islands, obtained 
from JSJ using the effective modulation amplitude (Afff). 

We also calculate the time-averaged longitudinal par- 
ticle distribution profiles in time domain, which is a his- 
togram of the particles position in time for a period much 
longer than the synchrotron one. This result simulates 
the profile that would be seen when measuring the signal 
from a button pickup, and is shown in Figure [5] 

Finally we simulate the damping of a coherent dipolar 
oscillation by phase modulation. The dipolar oscillation 
of a high current single bunch (30 mA) is excited by the 
LI HOM and, when phase modulation is turned on, the 
oscillation amplitude is reduced by a factor 30 (Figure 
I10II . which corresponds approximately to an attenuation 
of 30 dB. This is similar to the 40 dB damping of the dipo- 



lar spectral line produced by phase modulation (Figure 
|3J) that we typically observe in multibunch mode. Note 
that in this single bunch simulation the effects of an in- 
creased incoherent frequency spread are fully taken into 
account. 



C. Multibunch Simulations 

For the multibunch case, we simulate the dipolar os- 
cillation excited by the longitudinal HOM of the cavities 
fFigure lTlT \). The simulations are performed considering 
a total stored current of 200 mA. A comparison between 
two different multibunch tracking simulations shows that 
the change in the amplitude of the dipolar coherent syn- 
chrotron oscillation when the LI mode is turned on is 
greater than 100 dB, demonstrating the excitation of 
dipolar oscillations in the beam, as observed during users 
shifts. We could also observe a good agreement between 
the value of the growth time of the LI mode given by 
theory (4 ms) and simulation (4.5 ± 0.5 ms). 

Figures lllt \ and B show results of multibunch simu- 
lations without phase modulation for two different equi- 
librium bunch lengths. We see that even if the bunch 
length in increased by a factor of 3 (which is larger than 
the increase observed in multibunch simulations when 
phase modulation is applied) there is basically no damp- 
ing of the coherent oscillations excited by the LI HOM. 
In other words, the mere dilution of the beam bunches is 
not enough to explain the coherent damping efficiency of 
phase modulation: the internal rearrangement of longitu- 
dinal phase space and corresponding increased frequency 
spread is a fundamental ingredient, as could be seen in 
the single bunch simulations above. 

Finally, we simulate a multibunch case with phase 
modulation, which is shown in Figure ITTD . Now we ob- 
serve a damping of the CBM amplitude by a factor 5, 
which is much smaller than what is observed in the sin- 
gle bunch simulations as well as in experiments. This 
is a result of the hybrid character of our simulation, in 
which only one out of 148 bunches has an internal struc- 
ture (and corresponding frequency spread), although all 
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FIG. 11: Simulation of the centroid oscillation of a longitu- 
dinal CBM driven by an RF cavity HOM: (A) without phase 
modulation, (B) without phase modulation and three times 
the natural bunch length and (C) with phase modulation. The 
figure shows one of each 6 bunches, as it would be seen by a 
streak camera. 



bunches share the common increased bunch length. The 
damping of the structured bunch 60 weakens the coupling 
between the electron beam and the HOM reducing the 
amplitude of the CBM coherent oscillation. 



MEASUREMENTS AND RESULTS 



In Sections III and IV we calculate and simulate some 
characteristics of the longitudinal beam dynamics with 
RF phase modulation on the second harmonic of the syn- 
chrotron frequency. In this section we perform an exper- 
imental verification of the predictions made. The details 
of the experimental setup and the results obtained are 
shown in the following parts. 
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FIG. 12: General layout of the experimental setup. The solid 
line correspond to the experimental setup used to measure 
the bunch shape and the beam spectrum while the dashed 
line corresponds to the additional setup used to measure the 
damping time for coherent synchrotron oscillations. 



A. Experimental Setup 

The basic idea of the experiments is to excite the beam 
via phase modulation of the RF fields and observe its re- 
sponse as a function of a controlled parameter such as 
modulation amplitude and frequency. To modulate the 
phase of the accelerating fields, we use a voltage driven 
phase shifter attached to the output of the RF genera- 
tor responsible for producing the master RF signal. This 
phase shifter is controlled by a signal from another gener- 
ator that creates a modulation with a chosen frequency. 
The full experimental setup is shown in Figure IT2"1 

The shape of the bunch is measured by a wideband 
oscilloscope (HP 54750A, 20 GHz) and the amplitude 
of the dipolar synchrotron sideband is observed using a 
spectrum analyzer (HP8560) both receiving the signal 
from a button pickup in the ring. In the multibunch 
case we choose to observe the synchrotron line around 
harmonic ((281 of the revolution frequency that is close 
to the frequency of the LI longitudinal HOM from the 
cavity. 

We also measure the damping time for synchrotron 
oscillation of the centroid of a single bunch by com- 
paring the phase between the master RF signal to the 
476.066 MHz component of the longitudinal beam signal 
and observe the signal in a high sampling rate oscillo- 
scope (TDS684B). To create a longitudinal kick we use a 
function generator that produces a step function with a 
given amplitude and acts on a phase shifter in series with 
the one used to modulate the RF wave. The additional 
setup required to measure the damping time is shown by 
the dashed lines in Figure IT21 
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Single-bunch (5 mA) longitudinal profile observed in a oscilloscope with 
phase modulation (A^ = SO mi-ad) : 
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FIG. 14: Example of a measurement of the synchrotron damp- 
ing time. 



FIG. 13: Beam profile observed in an oscilloscope for a gap 
voltage of 460 kV. 



B. Experimental Results and Comparison with 
Theory 

In this section we present a comparison between the 
predictions of the theoretical model and measurements. 
This analysis is done in the following order: island for- 
mation as a function of modulation frequency, damping 
time of the centroid synchrotron motion as a function 
of modulation frequency and amplitude, measurement of 
island tune and its behavior as a function of modulation 
amplitude and frequency, and measurement of the BTF. 



1. Island Formation 

The time-averaged longitudinal bunch profile predicted 
by our single bunch simulation, Figure is compared 
with the experimental results of the signal from a button 
pickup observed on an oscilloscope and shown in Figure 
1131 We observe that, for modulation frequencies u> rn < 
u! s (2 — e/2), the profile corresponds to the case with three 
stable regions whereas, when u) m > uj s {2 — e/2), only two 
stable regions are created. This is an indirect way to 
observe the bunch internal structure created by phase 
modulation. The same result has already been obtained 
directly in other laboratories |ll| using a streak camera 
to observe the beamlets. 



Figure 1151 is the damping time ratio for a single bunch 
with phase modulation turned on and off. 

The results presented in FigurefHIare not the ones ob- 
served directly on the oscilloscope. When performing this 
measurement we observe that the phase of the RF compo- 
nent of the beam oscillates not only with the synchrotron 
frequency but also with multiples and sub-multiples of 
the modulation frequency. The raw data acquired is fil- 
tered and those frequencies removed from the spectrum 
in order to leave just a band of frequencies around the 
dipolar frequency of centroid oscillations, which allows 
us to determine the damping time of coherent dipolar 
oscillations. 

A set of measurements of damping time versus modu- 
lation frequency is shown in Figure 1151 and it is possible 
to observe that phase modulation is effective only for a 
narrow band of frequencies approximately 400 Hz wide. 
For frequencies above 50.2 kHz the effect no longer exists, 
which agrees with the theoretical model used to describe 
the damping mechanism. This model predicts that, for 
frequencies above 50.3 kHz (/„, = (2+e/2)/ s ), the islands 
collapse into one so that there is no increase in damping 
time. 

Using equation (|41|) , we can estimate the value for the 
damping time increase due to phase modulation and Fig- 
ure shows good agreement between theory and exper- 
imental data. The error bars are relatively large due to 
the difficulties of measuring the damping time. 



2. Longitudinal Damping Time 

The measurements of the damping of coherent syn- 
chrotron oscillations are performed using a function gen- 
erator and a phase shifter to apply a longitudinal kick 
to the beam. Using the signal from a button pick-up, we 
measure the phase of the RF component of the beam with 
respect to the master RF signal. Figure^jjis an example 
of the measurement for a single bunch with 4.3 mA and 



3. Measurement of the Island Tune 

An FFT of the damping time measurement data also 
shows a lower sideband of the line corresponding to half 
of the modulation frequency. This sideband, shown in 
Figure IPol correspond to the island tune [Appendix IB] . 
Figures IT7I and IT51 show the measured dependence of the 
island frequency with modulation parameters and the 
corresponding theoretical curves l|31f) . which are in good 
agreement. 
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FIG. 15: Increase in damping time as a function of modulation 
for a gap voltage of 460 kV: (red dots) measurements and 
(blue curve) theoretical predictions calculated for an island 
population of 50%. 
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FIG. 19: Experimental setup used for measuring the BTF. 
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FIG. 16: Observation of the island tune. 
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FIG. 17: Island frequency as a function of modulation fre- 
quency for a gap voltage of 460 kV. 



4- Measurement of the BTF 

The experimental setup used to measure the BTF is 
shown in Figure El We use a network analyzer to excite 
a dipolar oscillation on a low current single bunch and 
measure the beam response in amplitude and phase using 
a button pickup. 

The results for a set of measurements with phase mod- 
ulation turned on and off are shown in Figures l201andl?Tl 
respectively. The dots are the experimental data and the 
full line comes from the theoretical model. The result 
of those measurements indicates that the phase modu- 
lation produces an increase in frequency spread within 
the bunch since the width in the amplitude response it 
related to this spread. 

Figure 1221 shows the stability diagram for the same set 
of measurements above fFiguresl251and l2l"J> and the corre- 
sponding theoretical curves. The green dots and red line 
correspond to the case without modulation, whereas the 
black dots and blue line correspond to a phase modula- 
tion amplitude of 35 mrad and frequency of 51 kHz. Note 
that the stable area is effectively increased when phase 
modulation is turned on and also that some features, such 
as the loop predicted by the model, are present in the ex- 
perimental data. 
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FIG. 18: Island frequency as a function of modulation ampli- 
tude for a gap voltage of 460 kV. 



VI. DISCUSSION AND CONCLUSIONS 

We have shown experimental and theoretical evidence 
that, when RF phase modulation on the second harmonic 
of the synchrotron frequency is used, the available sta- 
ble area to the electron beam is effectively increased and 
that the mechanism responsible for this enhancement in 
stability is an increase in the incoherent frequency spread 
which, via Landau damping, helps to alleviate unstable 
CBMs. We also simulate and calculate the island forma- 
tion and perform measurements that indicate that these 
structures are indeed formed in the longitudinal phase 
space. 

Our simulations also indicate that the lengthening of 
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FIG. 20: Measurement of BTF for a single bunch without 
phase modulation for a gap voltage of 500 kV. The dots cor- 
respond to data points and the full red line to theory. 
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FIG. 21: Measurement of BTF for a single bunch, for a gap 
voltage of 500 kV, with phase modulation using the following 
parameters: f m = 51 kHz, A m = 35 mrad and N c =65 %. 
The dots correspond to data points and the full red line to 
theory. 
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FIG. 22: Stability diagram of a single bunch with (black dots) 
and without (green dots) phase modulation and the respective 
theoretical curves. 



diagram), although the island population still remains 
as a free parameter. We hope to be able to access this 
quantity in the future by using a streak camera. 

We have been using RF phase modulation since 2004 
to stabilize the electron beam of the LNLS storage ring. 
Although in the Brazilian machine the overall lifetime 
increase (we observed a lifetime increase of 30 %) was 
not as great as the one observed in other laboratories 
[!J [T(J, [HI the orbit fluctuation caused by the coupling of 
the beam with a HOM of the RF cavities was successfully 
suppressed improving the beam quality for users. 
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APPENDIX A: LONGITUDINAL HAMILTONIAN 
AROUND THE SECOND HARMONIC 
RESONANCE 



the bunches which results from phase modulation (an 
effect which has already been used at several laboratories 
as a tool to increase Touschek lifetime of intense bunches) 
cannot explain by itself the stabilization of CBMs that is 
observed when phase modulation is applied. In fact, only 
when bunches with internal structure (and corresponding 
frequency spread) are included in the simulations do we 
observe the intense damping of CBMs that is obtained in 
experiments. 

Finally, the model developed to describe the BTF mea- 
surements presents all interesting features of the exper- 
imental results, namely the double peaks and a second 
phase jump (which shows up as the loop of the stability 



Here we present the detailed derivation of the Hamil- 
tonian with phase modulation. Beginning with 0] and 
changing variables to 

4> = VzJ cos and 5 = — \Z2jsin-0, (Al) 
we get 

H(J,rp) — lu s J sin 2 -0 — uj s \/~2J tan(f> s costA 

+uj s tan cf) s sin[V / 2jcos ip + A m sin u> m t] 
— u> s cos[V / 2Jcos-(/' + A rn smLo m i\. (A2) 
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Expanding the terms in sin 'J' and cos 'J' using Bessel 
functions and restricting our analysis to even resonances 
close to the second harmonic of the synchrotron fre- 



J 



quency (u> m ~ 2tv s ), we can simplify the equation above 
to 



• cos 2ip ) J 
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■u! a tan tp^J 2Jcos tp + u s tan ip s V 2J sm(oj m t) Jq( V 2J) 



2w s ^(-l) fc J 2fe (V 2 J) cos(2fe^) + 2w s tan > s ^(-l) fc J 2fc+1 ( ^2 J)) cos[(2/c + 1) 



k=l 



fe=0 



-w s tan ^4A m J 2 ( v 2J) sin(o; m i — 2'0). 



(A3) 



r 



The new Hamiltonian is still time dependent and in or- 
der to eliminate this dependence we can perform a canon- 
ical transformation using the following generator function 



r 2i ./.,, = (v.-^-£).7 .A-h 



so that the new coordinates in a rotating reference frame 
are 



J = — - = J and tp = — — 

dtp oJ 



and the Hamiltonian is 

dF 2 



K = H(J, tp) 



dt 



H(J,f/>) 



(A5) 



(A6) 



performing a time average, the fast oscillation terms will 
vanish and we find that 



(Jf>« = ( W .-^)j- 



uj,J 2 uj,cJ 
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cos 2tp (A7) 



which is the Hamiltonian around the second harmonic 
synchrotron resonance. 



APPENDIX B: CALCULATION OF THE ISLAND 
FREQUENCY IN THE LABORATORY 
REFERENCE FRAME 

In the rotating reference frame the motion of the par- 
ticles in the islands can be described by 

<j>' = r cos(w'i) + xo and 5' — —r sin(a/t) (Bl) 

where tq is the amplitude of the oscillations around the 
stable fixed point, xq is th e sta ble fixed point coordinate 
in phase space and u>' = V AB is either the lateral island 
frequency or the central island frequency. To go from 
the rotating frame (which rotates clockwise with respect 
to the lab frame) to the laboratory reference frame, it is 



necessary to transform the above set of equations using 
the following rotation matrix 



sin(^) cos(^) 



n{uj m t/2) 

and the new particles coordinates are 

(*) 

or writing explicitly for each coordinate 

t 



To COS 



ro sin 



Xq COS 



Xq sin 



UJ m t 
UJ m t 



(B2) 

(B3) 

(B4) 
(B5) 



To understand this set of equations we have first to know 
wether particles in the islands rotate clockwise or coun- 
terclockwise. Taking the harmonic part of the Hamil- 
tonian (|21|l expanded around each stable fixed point we 
find that 



dd>' , ,, dS' 
— = AS' and — = 
dt dt 



(B6) 



so that if A > and B > the particles rotate clockwise 
and if A < and B < the particles rotate counterclock- 
wise. 

Calculating the coefficients A and B for our parame- 
ters, we find that particles in the central island rotates 
clockwise and particles in the lateral islands rotate coun- 
terclockwise. Using this result we can write the particles 
equations in the lab reference frame 

Ay 



(pi = r cos 
5 t = r sin 



Ui J t 

-Ui)t 



- Xo cos 
xq sin 



2 
2 



(B7) 
(B8) 



and 



So 



ro cos 
tq sin 



(^H< 

— +LU c )t 



(B9) 
(BIO) 
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for the lateral islands and central island respectively. 

From the equations above we note that when measur- 
ing in the laboratory reference frame the particles re- 
spond to and external excitation with three different fre- 
quencies: a narrow line at oj m /2, which is related to the 



to the external phase modulation, and a broad line at 
u> m /2 — u>i and uj m /2 + lo c that correspond to the parti- 
cles incoherent revolution frequencies around each stable 
fixed point. 
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